Azimuthal Anisotropy: Ridges, Recombination 
and Breaking of Quark Number Scaling 
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Azimuthal anisotropy is studied by taking into account the ridges created by semi-hard scattering, 
which is sensitive to the initial spatial configuration in non-central heavy-ion collisions. No rapid 
thermalization is required. Although hydrodynamics is not used in this study, the validity of hydro- 
dynamical expansion is not excluded at later time after equilibration is achieved. Phenomenological 
properties of the bulk and ridge behaviors are used as inputs to determine the elliptic flow of pion 
and proton at low px- At intermediate px the recombination of shower partons with thermal partons 
becomes more important. The (p dependence arises from the variation of the in-medium path length 
of the hard parton that generates the shower. The pr dependence of V2 is therefore very different 
at intermediate pt compared to that at low pr- Quark number scaling of V2 is shown to be only 
approximately valid at low pt, but is broken at intermediate pr, even though recombination is the 
mechanism of hadronization in all pr regions considered. 

PACS numbers: 25.75.-q, 24.85. +p 
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I. INTRODUCTION 

In relativistic heavy-ion collisions the subjects of 
transverse momentum {pr) distribution and azimuthal 
anisotropy have been given prominent attention both ex- 

Sirimentally and theoretically from the very beginning 
-0. In the second half of this decade details of the 
properties of elliptic flow have continued to be studied 
experimentally with greater accuracy but there 

is a waning of further theoretical development, at least 
in the light quark sector. The hydrodynamical model at 
low Pt [lJ]-pJZ| and the recombination-coalescence model 
at intermediate px [3- [Ml have described the data on 
elliptic flow so well that little room seems to exist for 
further improvement. In particular, quark number scal- 
ing (QNS) is a property that has been scrutinized to the 
point where it has been regarded as a strong evidence for 
the partonic degree of freedom before hadronization and 
for recombination at low px ■ In this paper we study the 
problem of azimuthal anisotropy in the framework of our 
version of the recombination model (2^ in which we con- 
sider thermal and shower partons. We take into account 
the ridges at low px generated by semi-hard scattering 
near the surface and the shower partons at higher px- 
We calculate the second harmonic V2 for pion and pro- 
ton, as well as for a light quark, and show that, when px 
is extended to the intermediate region, there is signih- 
cant departure from the QNS result suggested by simple 
consideration of quark recombination [20l . [2l| . 

We state at the outset that we do not have a dynam- 
ical description of the time evolution of the expanding 
system. We do not use hydrodynamics explicitly, but 
by considering thermal distribution of parton just be- 
fore hadronization we imply the validity of hydrodynam- 
ical expansion at some point in the evolutionary history. 



However, we do not subscribe to the applicability of hy- 
drodynamics at early time, such as at tq = 0.6 fm/c 
p^ . [l7| . Rapid thermalization has not been shown to be 
the consequence of any dynamical process that is firmly 
grounded and commonly accepted. 

It has been pointed out in Ref. [1^ that there exists 
an alternative mechanism to relate the spatial asymmetry 
at early time to the momentum anisotropy at late time, 
without relying on the assumption of fast equilibration. 
That mechanism is semi-hard scattering, soft enough to 
have high probability of occurrence, but hard enough to 
take place at r < 0.2 fm/c. The phenomenology asso- 
ciated with such processes may be termed "ridgeology," 
which is the study of the properties of ridges that have 
been found to accompany jets, even when the initiating 
jets are weak and the peak-to-ridge ratio is small [24, 25,]. 
In this paper we develop further the study of elliptic flow 
based on ridges at low px and shower partons at inter- 
mediate Px- 

The shower partons used in our model are defined at 
the hadronization scale so that the recombination of two 
types of shower partons in the same jet reproduces the 
fragmentation functions (FF) of a hard parton to a spe- 
cific meson [26]. Similarly, three quarks in the shower 
can recombine to form a baryon in agreement with the 
baryon FF without adjusting any free parameters [27| . 
A shower parton can also recombine with a thermal par- 
ton in the immediate vicinity of the jet (but not belong- 
ing to a part of the jet), both being soft enough to un- 
dergo hadronization. It has been shown in [22] that such 
thermal-shower (TS) recombination predominate in the 
region 3 < px < 6 GeV/c. Our concern for V2{px) in this 
paper will be for < pt < 6 GeV/c. 

Above Px = 6 GeV/c, SS (or SSS) recombination be- 
comes important; that process is equivalent to fragmen- 
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tation, since that is how the shower parton distributions 
are determined in the first place [26]. With that duahty 
of fragmentation and recombination in mind, it is easy 
to see that QNS cannot be vaUd at high px, since jet 
quenching responsible for V2{pt) involves only one hard 
parton, not two (or three) semi- hard partons. Although 
we do not consider the region > 6 GeV/c, there is still 
the question: at what point does QNS begin to break 
down. Our study here shows that the breaking of QNS 
begins at the transition point between TT and TS re- 
combination for pion and between TTT and TTS-fTSS 
for proton. The reason is in the nature of TS and TTS 
recombination, a bad approximation of which can lead 
to a simplistic formula that falsely suggests QNS. Data 
that seem to support QNS are all for minimum bias, and 
do not reach the upper region of intermediate px- 

In Sec. II we give the general formulation of how az- 
imuthal anisotropy can arise from ridges without any de- 
tails on the pt dependence. Elliptic flow is then calcu- 
lated for pion and proton production at lowpr in Sec. Ill, 
where only thermal partons are considered. In Sec. IV 
the study is extended to intermediate px-, where the con- 
tribution from shower partons is included. The breaking 
of quark number scaling is investigated in Sec. V. The 
final section contains the conclusion of this work. 



II. AZIMUTHAL ANISOTROPY ARISING 
FROM RIDGES 

Let us first review our approach to elliptic flow at low 
Pt without rapid thermalization 23]. Instead of assum- 
ing the meaningfulness of thermodynamical quantities, 
like pressure and temperature, at early time, we recognize 
that hard scattering of partons can occur at all virtual- 
ity Q^, with increasing probability at lower and lower Q^, 
and that when the parton transverse- momentum, kx, is 
around 2-3 GeV/c, the rate of such semi-hard scattering 
can be high, while the time scale involved is low enough 
0.1 fm/c) to be sensitive to the initial spatial con- 
figuration of the collision system. When such scattering 
occurs near the surface of the overlap region in the trans- 
verse plane, each semi-hard parton creates a ridge in Arj 
and A0 |29|]. When triggers are used, ridges have been 
found in the associated-particle distribution on the near 
side with trigger momentum p^'^ > 4 GeV/c [2^. More 
recently, p^'® has been reduced to as low as 2.2 GeV/c 
and p^'^°'^ as low as 1.5 GeV/c, where the ridge yield 
significantly dominates over the yield of the peak that 
sits above the ridge [1^. It suggests that ridges due to 
the scattering of low-x partons (< 0.03) are abundantly 
produced even if triggers are not used to select events 
to examine their properties. The effect of such ridges on 
both the px and (f) dependences of the produced particles 
should not be ignored. 

If the semi-hard scattering occurs in the interior of 
the dense medium, the energy of the scattered partons is 
dissipated in the medium and contributes to the thermal- 



ization of the bulk. That process may take some time to 
complete, a likelihood that is acceptable in the approach 
adopted here, since we have no need to require thermal- 
ization to be fast. If the semi- hard scattering occurs near 
the surface of the medium, one of the scattered partons 
may be directed outward and lose energy to the medium 
on its way out. The enhanced thermal partons near the 
jet trajectory can recombine and form hadrons in the 
ridge. This interpretation of the ridge has been applied 
successfully to triggered events [s^], and provides a res- 
olution to the O puzzle [3l|, [13] ■ The same mechanism 
of ridge formation is, however, also valid, if no trigger is 
used. Although the variables A77 and Acj), usually defined 
relative to the trigger momentum, would be meaningless 
without a trigger, the existence of ridge in the difference 
variables t^a and 4>a in autocorrelation that used no trig- 
gers has been well established [2^, . 

The direction of a scattered parton is random, but af- 
ter it is averaged over all events, the average direction 
of all outward partons near the surface is normal to the 
surface. Since low-x partons are copious and the rate 
of semi- hard scattering with fcy < 3 GeV/c is not sig- 
nificantly suppressed, the semi-hard emitters can form a 
layer along the surface of the overlap region prescribed 
by the geometry at the initial time of collision. The for- 
mation of the ridge of hadrons takes some time for the 
transverse expansion to complete, but the directions in 
which the ridge partons flow are determined by spatial 
configuration at early time. We do not rule out the appli- 
cability of hydrodynamics at some point of the expansion 
process when equilibration is established. However, fast 
thermalization is not needed if semi-hard scattering can 
initiate the anisotropic expansion. 

Details about ridge formation are rather complicated, 
especially if a model is to be successful in reproducing the 
very recent data on the ridge yield as a function of the 
trigger azimuthal angle relative to the reaction plane ^3] . 
A full treatment of that problem is a separate subject of 
its own [ssj . and is unsuitable for inclusion here; further- 
more, we do not use triggers in the study of inclusive dis- 
tributions from which V2 is to be determined. However, 
aspects of that problem are needed to demonstrate the 
process of averaging over the jet direction of the semi- 
hard parton. We describe in the Appendix our model 
calculation of the ridge distribution in the azimuthal an- 
gle (j), and show how the result can be represented by a 
simple box approximation. We take that approximation 
as the starting point here and proceed to the calculation 

of V2{px)- 

The overlap region in the transverse plane for two nu- 
clei of radius Ra at impact-parameter h apart is, assum- 
ing simple geometry with sharp boundaries, the almond- 
shaped area bounded by two circular arcs whose maxi- 
mum angle is $, where 

cos$ = 6 = 6/2i?A , (1) 
and the angle within the arcs satisfies 4> £ TZ, which is 
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a set of angles defined by 



101 < $ and Itt - 01 < $. 



(2) 



It should be emphasized that the almond-shaped region 
is relevant for the consideration of the initial problem of 
semi- hard scattering whose time scale is short, although 
hadrons in the ridge are formed later when the elliptic 
geometry is more pertinent. In the Appendix we show 
that the ridge (R) distribution in <p can be represented 
by the box approximation 



where 



6(0) =6l($- |0|)-F6l($- |7r-0|) 



(3) 



(4) 



The bulk (B) medium has no dependence and will be 
denoted by B(jjt)- The single-particle distribution at low 
Pt is then 



dN 



PTdprdip 



= B{pt) + R{ptMcI^) 



(5) 



This is a simple expression in closed form that enables 
us to calculate V2{pt) analytically. 

The second harmonic in the distribution is 



V2{pt) = (cos 20) 



Jp ^ d<j) cos 20 dN/prdpj 



Jg ^ d(j) dN/prdpTdc/) 
When Eq. ^ is used in the above, we obtain 



V2{PT,b) 



sin2$(6) 



ttB{pt)/R{pt) + 2$(&) 



(6) 



(7) 



If the first term in the denominator is much larger than 
the second, as we shall show below for low pt-, then we 
have the even simpler formula 



V2{pT,b) 



R{vt) 

■kB{pt) 



sin2$(&) 



(8) 



where the px and h dependences are factorized. Thus the 
centrality dependence at fixed pt is essentially specified 
by sin2$(5), and is in accord with the data at low pT 
[23j . We shall show more detailed properties of the cen- 
trality dependence below. Equation ([8]) is the analytic 
result that represents in a simple approximation the con- 
sequence of considering the effects of semi-hard scattering 
instead of fast thermalization. 



III. ELLIPTIC FLOW AT LOW pr 



at low qT- This distribution includes the effect of semi- 
hard partons created in the interior of the medium; they 
lack the energy to reach the surface to get out. The en- 
ergy lost to the medium is thermalized and contributes 
to the value of the effective temperature T. The time it 
takes for the weak jets to thermalize need not be short; 
it can take a significant part of the expansion phase of 
the whole system, if necessary. If the semi-hard scatter- 
ing occurs near the surface, one of the scattered partons 
may emerge, while the recoil parton directed inward gets 
thermalized. The emerging semi-hard partons along the 
surface lose extra energy to the medium in addition to 
those others that cannot escape. Thus there is an en- 
hancement over the bulk, for which the thermal distribu- 
tion has the same form but with a higher inverse slope 

r 



dgrdq) 



e n 



(10) 



The difference between Eqs. ([9]) and (jTO]) is the ridge ef- 
fect. Note that the normalization factor C is the same 
in Eqs. © and (fTUl) because there is no difference be- 
tween the weak jets that fail to emerge from the surface 
and those that barely emerge with negligible energy to 
develop any additional enhancement at qt = 0. It is only 
when the semi-hard scattering occurs sufficiently near the 
surface that a substantial ridge can be formed with non- 
vanishing qx- The two equations above are not derived 
from any dynamical equation of time evolution, but are 
to be inferred from the data on particle distribution in 
Pt after the quarks hadronize by recombination. Clearly, 
the use of thermal distributions implies equilibration be- 
fore hadronization, but it need not be accomplished at 
early time. It is not necessary for us to specified how 
long the equilibration time is, since Eqs. ([9]) and (fTO]) 
may be regarded as phenomenological input with T and 
T' to be determined from data. 

To derive the pion and proton distributions from Eqs. 
@ and (llOp the formalism in (2^ is to be used. How- 
ever, in [22] only central collision is considered, for which 
the quark distributions are assumed factorizable before 
recombination. Here for the study of anisotropy all 
centralities must be considered, and the assumption of 
factorizability of uud distribution for the production of 
proton in peripheral collision (where thermal parton den- 
sity is low) is questionable. We shall treat such compli- 
cations in a later section. For now, let us ignore the issue 
of non-factorizability and proceed as in (22j for central 
and mid-central collisions so as to show the connection 
between the ridge effect and elliptic flow. 



Let us now consider the low-px behaviors of B{pt) and 
R{pt)- The parton distribution in transverse momentum 
qT of the bulk medium has the thermal behavior [22l | 



Pion 



(9) 



Starting from Eq. using valon distribution in the 
recombination function, and neglecting pion mass, we ob- 
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tain in |22l | the pion distribution due to TT recombination 



B^{pt) 



dNi 



prdqj 



6 



'Pt/T 



(11) 



for the bulk medium at any (j). The normahzation factor 
will be canceled at low pr, so it need not be specified 
here. T will be discussed below. It should be recognized 
that the form in Eq. ^ is chosen so as to yield the ex- 
ponential distribution in (jlip . Starting from Eq. (IIOD we 
obtain similarly for (j) G TZ 



B-^{pt) + Rtv{pt, 



prdqrdcj) 



6 



-Pt/T' 



(12) 



While Eqs. pT|) and are consequences of ^ and ([TT 
for any pt^, their phenomenological application is useful 
only for low px where recombination with shower par- 
tons is relatively unimportant. However, for pT ~ Q the 
hadronic mass cannot be ignored, even for pion, let alone 
proton. The recombination model we use is based on a 
formalism for large momenta [36| - [38| . where momen- 
tum fractions are meaningful. Thus it is an assumption 
that the model remains quantitatively valid at low pT, 
an extrapolation that can be made more acceptable if 
the mass effect can be taken into account. To that end, 
we adopt the ansatz that pT in Eqs. (fTTjl and (fT2)) are to 
be replaced by the transverse kinetic energy Et 



Et{pi 



rriT - mo, 



rriT 



2\l/2 



where toq is the hadron rest mass. We then have 

C72 



-Et[pt)/T 



6 



(13) 



(14) 







-Et(pt)/T' 



, <t>en, (15) 



where T and T' remain to be determined from data. Sub- 
tracting Eq. from pS]) results in ([3]), where R.^{pt) 
has now the concrete form 

i?^(p^)^^e-^-(P-)/^'(l-e-^-(f-)/^") , (16) 



11 1 _ AT 



AT = T' -T 



(17) 



For the values of T and T', we examine the data from 
STAR that has focused on ridgeology [11]. It has been 
found that the ridge distribution in p?^^""^ is nearly ex- 
ponential and remains essentially the same for a wide 
range of trigger momentum p^'® > 4 GeV/c. The min- 
imum value of p?^^°'^ measured is, however, 2.2 GeV/c, 
too high to exhibit the small Et behavior in Eq. 
Furthermore, it is also shown in Ref. that the p/tt 
ratio in the ridge for 4 < p^'^ < 5 GeV/c is about 4.5. 



It means then that the data on ridge in [24| for unidenti- 
fied charged particles cannot be used to determine T' in 
Eg. (|T6l) for pion. It should also be pointed out that in 
[24| the ridge slope is compared to the inclusive slope of 
charged particle distribution for the same range of pt as 
for 2.2 < pf^'^°'' < 4 GeV/c. That is the region in which 
the p/tt ratio of the inclusive single-particle distributions 
is approximately 1. So no information about the pion T 
for the bulk can be obtained from Ref. [24[ either. 

For central collisions we have 6 ~ and $ ~ 7r/2, 
so B{pt) + R{pT,(p) is for all (p without anisotropy. 
The single-particle distribution is therefore a measure of 
B + R, not B alone. Semi-hard scattering is present in 
all events so ridges contribute to dN/prdpT whether or 
not triggers are used to select jet events. Thus Eq. fTS]) 
can be applied to the inclusive data for identified pions 
to determine T' . The data are not exactly exponential 
for all low Pt- For < 1 GeV/c there are pions from 
resonance decays. We extract the value of T' from the 
region px > I GeV/c, in which the ridge particles from 
the enhanced thermal source make definitive contribu- 
tion. Using the data in Ref. [1] we fit the 7r+ distribu- 
tion in Et at 0-5% centrality for 1 < Et < 3 GeV and 
get T' = 0.3 GeV. If TT recombination were dominant 
in peripheral collisions, we need only go to the tt"*" distri- 
bution at 80-92% centrality and determine T, since the 
ridge contribution is small when <i> is small. But that 
is not the case. At large b the thermal source is weak, 
and SS recombination dominates. It is shown in Ref. (39j 
that SS > TT for pT > 2 GeV/c at 80-92%, as in pp col- 
lisions. Without direct experimental guide to determine 
T for Bt^{pt), we proceed by adopting AT = 45 MeV 
in common with what we shall be able to obtain in the 
study of the proton case to be described below; it is con- 
sistent with the ran ge o f values suggested experimentally 
in [Hi and used in j23l|. Using Eq. p7)) we then obtain 
T ^ 0.255 GeV and 



T" = 1.7 GeV 



(18) 



for the pion. 

The ratio B.^{j)t) / Rtt{pt) has a simple form when Eq. 
is rewritten as 

i?^(pj,) ^ ^e-^-(P-)/^ (e^-(f-)/^" - l) , (19) 
so that the prefactors cancel and we have 

(20) 



1 



B.{pt) ^ 



which depends only on T" . We assume the validity of this 
equation for all pT <2 GeV/c where TT recombination 
is valid. When pT is small, Eq. can be approximated 
by T" /Et, so Eq. ([8]) has the simple expression 



^2bT,6) = ;^Msin2$(6). 



(21) 
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FIG. 1: (Color online) Comparison of calculated with data 
for Au+Au collisions at 200GeV 0] for 8 centrality bins whose 
corresponding values of b are shown in the legend. 



The initial slope is 



Et=0 



ttTI 



■sin2$(fe) 



(22) 



which is dependent on only T" apart from the geometri- 
cal factor sin2$Y&). That property is in good agreement 
with the data [3], as shown in Fig. 1. Equation ((22)) is 
a distinctive feature of V2 that is driven by ridges. For 
non-vanishing values of Et we use the full expression of 
Eq. (PO)) in ([7]), obtaining the result shown in Fig. 1 for 
Pt < I GeV/c. Deviation from Hnearity is perceptible. 
There is some discrepancy between our results and the 
data, the most noticeable region being around E^ ~ 0.65 
GeV, but taken as a whole the overall agreement with the 
data is good. That is remarkable, since we have not ad- 
justed any free parameters in order to fit. 

It is worthwhile to digress and explain why our model 
that depends on the initial geometry of the medium can 
produce satisfactory V2{pt), whereas earlier model based 
on surface emission without hydrodynamical expansion 
failed [40]. In our case, although the average direction of 
the semi-hard partons are emitted normal to the surface 
of the initial configuration, the hadrons in the ridges are 
produced after expansion. The late emission of hadrons 
introduces a scale in Rtt{pt) / B.^{pt) in the form of T", 
which controls the normalization in Eq. ((22l) . The early 
emission of semi-hard partons determines the (f> depen- 
dence through (E>(6) in agreement with data. In the case of 
po| the only scale is the temperature T = 140 MeV with- 
out expansion, so V2{pt) saturates at around pT ~ 0.2 
GeV/c in disagreement with data. Furthermore, it is 
not detailed enough to give centrality dependence. The 
importance of hydrodynamical expansion is pointed out 
in [40| . on which we have no disagreement. In sum- 
mary, ridge physics is not a simple surface-emission prob- 
lem because it comprises both early-time parton emis- 
sion and late-time hadron formation, and therefore has 



more dynamical content to result in the proper azimuthal 
anisotropy. 



B. Proton 

For proton production in central Au-Au collision we 
have obtained for TTT recombinations 1231 



dNp 
Prdpi 



Po 



■Pt/T 



where 



A 



B{a + 2, 7 2)B{a + 2, a + 7 4) 
6 B(a-f l,7-Kl)B(a + l,a + 7-l-2) 



(23) 



(24) 



a = 1.75 and 7 — 1.05. These beta functions come from 
the recombination function of proton, which depend on 
valon distribution characterized by a and 7 [411] . The 
factor p^ arises from the integration over the momenta 
of the quarks that recombine. The factor pq comes from 
the invariant distribution p^dN/d^p on the left-hand side. 
When applied to pr > 2 GeV/c at mid-rapidity, the 
Pt/po fs-ctor was approximated by pr in ^22]. Here we 
want to extend Eq. ([23]) to lower pT, still at y « 0, so to 
take the mass effect into account we rewrite the equation 
in the form 



niT 

Similarly, for bulk + ridge we have 



(25) 



BpipT) + RpiPT,^) = A^e-'^-^P-^/^', ^en. (26) 
rriT 



The ridge solution is then ior (f> £ TZ 



Rp{pT) = A^e^ 
rriT 



Et(pt)/T' 



(27) 



When the exponential factor is in terms of pr, as shown 
in Eq. ((23l) , the value of T should be the same as that in 
(fTTj) for pions, as well as that in © for the quarks. That 
is because of the 5 (pt — I^t) function in the recombi- 
nation function that preserves the exponential behavior 
from quarks to pions to proton. Experimentally, it has 
been found that the pr distributions of tt"*", if"*" and p 
have the same inverse slope in their exponential behavior 
for 1.5 < Pt < 3 GeV/c (see Fig. 5 of [^|), thus confirm- 
ing the result of the recombination model. However, the 
same data, when plotted in terms of Et, show different 
slopes (see Fig. 9 of [1]). As stated in the Sec. HI. A 
above, the slopes determined from the Et distribution 
for 0-5% centrality is the value of T'. We obtain from 
the proton distribution in Q = 0.35 GeV. It is now 
possible to extract some information from the ridge dis- 
tribution in Ref. which is for 0-10% centrality. The 
lowest p!^'® range is 4 < p^'^ < 5 GeV/c for which the 
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p/tt ratio is 4.5. Assuming dominance by proton, we use 
Eq. (P7|) to fit the data by varying AT. Since the nor- 
malization of the ridge distribution in [23] is dependent 
on p^'®, whereas Eq. (P7|) makes no explicit reference to 

p^'®, we focus only on the pt dependence of the par- 
ticles in the ridge and adjust the normalization to fit. 
The best fit is for AT = 45 MeV, for which the result is 
shown in Fig. 2. The data points are outside the pr range 
where the dip occurs at low pT\ nevertheless, a good fit 
is achieved in the region where data exist. It should be 
recognized that the factor p^ in front of exp{~ET /T') in 
Eq. (1211), as well as the factor, 1 - e:iq){-ET /T"), after 
it together make the effective inverse slope to be greater 
than T', where the data points are. A measurement of the 
dip in the proton distribution in the ridge at small pt in 
Fig. 2 will serve to test the validity of our model. With 
AT = 0.045 GeV we obtain from Eq. (HZl) T = 0.305 
GeV, and 



rpll 
P 



2.37 GeV 



(28) 



for 0-10% centrality. 



10" 




.4<p:^'y<5 GeV/c, 0-10% 
-R (p^) AT=45 MeV 



2 3 
Pt (GeV/c) 



FIG. 2: Fit of the pr distribution of charged particles in the 
ridge 24] by Eq. (p7l) . 

For centrality dependence we go to Ref. [1] and find 
that, whereas the slope for pion is essentially independent 
of c (defined as centrality in %), that for proton (and 
antiproton) decreases with c. Since the Et distribution 
of proton has a break in slope from Et < 1 GeV to > 1 
GeV, we choose to consider the tabulated slope for p, 
which is independent of the Et regions and in our view 
should be the same as for p. We find that (identifying 
T' = T'-) 

= 0.35 (1 - 0.5c) GeV, c = c/100 (29) 

gives a good fit of the data as shown in Fig. 3. The 
reason for the centrality dependence is that for baryon 
production the three quarks (or antiquarks) needed for 
recombination become less independent of one another at 



larger c, where less thermal partons are generated in the 
nuclear collision. We know that in very peripheral colli- 
sion the system should approach that of pp collision and 
few thermal partons are created for baryon production. 
When the uud joint distribution is not factorizablc, our 
model calculation of p production that is based on fac- 
torizablc quark distribution cannot be applied directly. 
Equation (|29p is the simplest way to summarize the c 
dependence of the inclusive distribution. 
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FIG. 3: Fit of inverse slopes of p production for five centrality 
bins H. 



Fixing AT at 45 MeV, we have 



Tp = 0.305 



0.175c GeV 



(30) 



Tp follows simply from Eq. ^T7\ . The result can be well 



approximated by a convenient formula 



rpll 

p 



2.37(1 - 1.05c -K0.26c2) GeV 



(31) 



The consequence on V2 can now be calculated as in the 
case of pion, using Eq. for low Et- The initial slope 
for proton is then 



DEt 



1 



Et=0 



7rT;(c) 



sin2$(&) , 



(32) 



Since T^' ^ T^', the initial slopes for proton and pion 
are not the same. That seems to violate quark number 
scaling, since dividing V2 and Et by the same does 
not change the initial slopes. However, the empirical ev- 
idence for QNS is only for minimum bias data [1, |^ [l2| • 
If we evaluate T^'(c) at c = 30%, we find T^'(0.3) to be 
very nearly 1.7 GeV in agreement with Eq. fTSl) . Thus 
our result is in accord with the experimental evidence for 
QNS at low Et, but also indicates violation of QNS in 
central and peripheral collisions. At higher Et there are 
other reasons for more serious breaking of QNS; that will 
be discussed in a later section. 

At higher Et, but still < 1 GcV, B^{pt)/ R^{pt) has 
the same form as Eq. ^0^ . so V2{pT,b) can be calcu- 
lated using Eq. 0. The result is shown in Fig. 4 in 
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FIG. 4: (Color online) Comparison of calculated iij with data 
for Au+Au collisions at 200GeV 0) for 8 centrality bins whose 
corresponding values of b are shown in the legend. 



rough agreement with the data that have large errors . 
Note that whereas (pt, b) saturates at around 40-60% 
centrality, V2{pT,b) continues to rise at more peripheral 
collisions. That is because, as b increases, Tp'(c) mono- 
tonically decreases so that sin 2$ (6) is not the only fac- 
tor that determines the dependence on centrality, as is 
the case with pions. The maximum of sin2<i>(&) occurs 
at $ = 7r/4, corresponding to 6 = l/-\/2, or c = 1/2. 
Equations (O, (pOj) and ([5^ indicate that the decrease 
of Tp (c) causes v^iEr, b) to continue to increase beyond 

b — 0.7 until the severe decrease of sin2$(6) at large 
h brings it down. We show in Fig. 5 the b dependence 
(for 2Ra = 14.7 fm) of v^{pT,b) and v^{pT,b) at fixed 
Pt = 0.52 GeV/c. The characteristics of the data Q are 
well captured by the simple formula, Eq. ([7]), for both 
pion and proton. 
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FIG. 5: Dependence of V2{pT,b) onb = b/2RA at pr = 0.52 
GeV/c for pion and proton. Data are from Ref. (3]. 



IV. AZIMUTHAL ANISOTROPY AT 
INTERMEDIATE pr 

As Pt is increased to above 2 GeV/c, it is neces sary 
to consider the role played by the shower partons [23 |. 
Both the ridge and shower partons are effects of jets, the 
former due to semi-hard partons, and the latter due to 
harder partons. The transition from one to another is, 
of course, a continuous one. Our recombination model 
was formulated in [2^ at a time before the ridges were 
discovered. In that model the three types of recombina- 
tion are TT, TS and SS, where no ridges are considered 
in the thermal partons. We now realize, as discussed in 
the preceding section, that the thermal distribution is 
B + R, since semi-hard scattering is always present. The 
only effect of this realization is just to relabel T in previ- 
ous work by T' now, as its value is determined from data. 
As we proceed to consider TS recombination in this sec- 
tion, it is T' that we shall use for the thermal partons. 
The condition of </> e 7^ in Eq. ^ for T' to be used is 
for TT recombination. Now for TS recombination, there 
is a hard parton to generate the shower parton. That 
hard parton may have any </>. Even if ^ 7?., the en- 
ergy loss of the hard parton can enhance the thermal 
partons near its trajectory, so the new T' can depend 
on b and (f>, characterizing the new ridge associated with 
hard parton. The dependence of the thermal partons 
is different from that of the shower partons. The latter is 
due to jet quenching of the hard parton that depends on 
path length in the medium. This change of the origin of 
the (j) anisotropy, when TT dominance is replaced by TS 
dominance at higher pT, is the basic cause for the change 
of the nature of pt dependence of V2{pt)- 

It should be noted that in other recombination models 
the azimuthal anisotropy at intermediate and higher px 
has been studied already [H, [l^ . Our approach here is 
different in the use of shower partons and in the inclusion 
of medium effect by treating TS recombination, whereas 
earlier studies considered fragmentation of hard partons 
as an additive component. The effect of jet quenching 
alone on V2 for pT > 2 GeV/c in a simple geometrical 
study of the (j) dependence is known to be too low |4J] . 



A. Hard Parton's (f> Dependence 

The shower parton distribution Sj (z) is the invariant 
probability of finding a parton of type j with momentum 
fraction z in a shower initiated by a hard parton of type 
i. Its detail properties are described in Refs. [l^jH^I, and 
its application to TS recombination in central collisions 
averaged over all </> is discussed in Ref. . We now con- 
sider (j) dependence due to energy loss of the hard parton 
with varying path length in the dense medium. Let us 
denote the distribution of hard parton i emerging from 
the surface of the medium with transverse momentum k 
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at angle by 



kdkdydcj} 



(33) 



Jet quenching degrades the hard-scattering momentum 
from the value k' at the point of scattering to the emerg- 
ing momentum k by an amount Afc that depends on the 
path length in the medium. The maximum length, 
^max of a straight line at angle </> that passes through the 
origin of the almond-shaped overlap region in the trans- 
verse plane satisfies the equation 



[In 



. cos 



+ 6) + 



■ sm( 



(34) 



Thus with the definition I — ^,„ax/2-RA we have 

i(6,0) = -6COS0+ (1 - S^sin^ 0)1/2 _ (35) 

This is a normalized path length that quantifies the de- 
pendence on b and <j). Assuming that the energy loss is 
proportional to the square root of the initiating parton 
momentum llOl l43ll. we write Afc in the form 



Afc = e{b)i{b, 0)a 



(36) 



where e(6) is the e nerg y-loss coefficient that may be given 
a reasonable form [131 



<h) = eo 



1 - 6-2(1-'') 



(37) 



since the density decreases with increasing b. The coeffi- 
cient eo is to be determined below. 

In Ref. [1^ we have found that the shower parton dis- 
tribution at mid-rapidity in a central heavy-ion collision 
is given by 



SAq)^(Y. I dkkMk)Si{q/k) 



(38) 



where /i(fc) is the distribution of hard partons i without 
nuclear suppression and ^ = 0.07 is the suppression fac- 
tor that is necessary to fit the normalization of the pion 
spectrum at intermediate px by TS recombination and is 
in accord with Raa being ~ 0.2 at higher px- We must 
now generalize Eq. ([55)1 to non-central collision. When 
6 > 0, there is dependence in Afc, given by Eq. ([5S[) . 
Instead of the suppression factor ^, we can describe the 
energy loss of the hard parton by shifting the parton mo- 
mentum and writing Sj (q) as 



J dkkmSlQ (39) 

i 

= J2 I dkkMk + Ak)Si (^^) 



at 6 = 0, and then generalize the RHS to 6 > and 
endowing it with a dependence through Eqs. l[55[) and 



([36[) . We note that fc is the momentum of the hard parton 
that emerges from the medium, and q is the momentum 
of the parton j in the shower, so fc is to be integrated 
from the lower limit at q. The momentum k' — k + Afc 
is that of parton i immediately after hard scattering and 
before traversing the medium. 

Writing Eq. ([55[) as fc' = fc -|- ei^/k', we can solve for 
k'{k,b,(f>), which, when substituted into fi{k'), yields a 
function of fc that is denoted by Fi{k,b,(j)) in Eq. (I33p . 



/,(fc'(fc,5,0)) = 2^F,(fc,6, 



(40) 



where Fi is defined per radian, while fi{k') is integrated 
over all 4>. The distribution fi{k') has a power-law be- 
havior 



(l + k'/kof 



(41) 



where the parameter a, fco and f3 are given in Ref. (4j| 
for a variety of parton type i. Thus Fi{k,b,4>) can be 
written in the form 



F,(fc, 6, 0) = ^/,(fc)G(fc, &,(/.) , 



where 



G(fc, &,(/.) = 



Ak{k,b,q 
fc + fco 



(42) 



(43) 



Keeping only the zeroth and second harmonics of the 
dependence, we have 



G(fc, 6, (j>) = go{k, b) + 2g2(fc, b) cos 20, 



(44) 



where go and 52 can be determined explicitly in terms 
of eo and b, beside fc, as we shall show in subsection C 
below. 



B. TS-I-SS Recombination 

Having obtained the (j) dependence of the hard parton 
distribution per radian at the medium surface, Fi{k, b, 0), 
we can proceed to the shower partons for recombination 
with the thermal partons. We have 

5,(gi,6,0) =^ J dkkF,{k,b,4>)Si{q^lk) , (45) 



7^.(q2,^<^) = C(6)g2e-«^/^'(^'^) , (46) 



]^ rPT 

— dq2 T{q2,b,(f))S{pT - q2,b,(j>),{A7) 
Pt Jo 



prdpTdcj) Px 

where the quark types j and j' are paired in appropriate 
ways to form a pion of a specific charge; for example for 
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7r+, they are ud and du. The pion recombination function 
has already been taken into account in the derivation 
of Eq. gl]) [23. The values of C{b) are given in Ref. 
[3^. The thermal parton distribution T{q,b,4>) is i3 + 
R with a new inverse slope T'[b,4') that depends on b 
and (j) because it is a measure of the enhanced thermal 
partons in response to the hard parton with momentum 
k and angle (f). Equation (|T7)) should not be considered 
for pt < 2 GeV/c, since the validity of the formalism 
is questionable at low and TS is dominated by TT 
recombination anyway. 

Unlike the (j) dependence at low px where semi-hard 
partons are restricted to €! 7?,, now the hard parton 
can originate from the interior of the medium and can be 
directed at any (j). The thermal partons in the vicinity 
of the trajectory are enhanced in proportion to the local 
medium density which is dependent on 6 and cj). Since 
that density at fixed < 6 < 1 decreases with increasing 
0, we adopt the reasonable form 



f '(6, 0) = T' 



\ + ab cos (j) 
1 + 2ab/TT 



(48) 



whose average over O<0<7r/2is the value T' = 0.3 
GeV used in Sec. III. A. The size of the parameter a awaits 
revelation by data on the ridge dependence on the trig- 
ger azimuthal angle, about which there is preliminary 
supportive indication from STAR. We shall use the pro- 
visional value a = 0.1 to carry out the calculation in Sec. 
IV.D. 

For SS recombination it is important to recognize that 
only one hard parton is involved and that the two shower 
partons are from the same jet. Thus the formalism is 
the same as described in Ref. [13], the only difference 
being the replacement of £,fi{k) /27r by Fi{k, b, (f>) for non- 
central collision, i.e., 



djVf (b) 
PTdpTdcj) 



— / dkkFi{k,b,(i 

X / dq{S,S}{q,k,pT) 
Jo 



(49) 



where the curly brackets denote the symmetrization of 
the leading parton momentum fractions zi =^ q/k and 
Z2 = (pt - q)/k 



{S, S}{q,k,pT) = ^ 



S: 



Sl{zi)Sl 
zi 



Z2 



1-zi 



1 



Z2 



Si {Z2) 



(50) 



As explained in Ref. [2^, the last integral in Eq. (|49|l is 
essentially the fragmentation function D{pT/k), so the 
equation describes both recombination and fragmenta- 
tion, neither of which has (/) dependence. Clearly, QNS 
cannot be valid when the only (p dependence arises from 
the hard parton distribution Fi(k, 6, 0). Let us write the 



sum of Eqs. p?)) and (1491) symbolically as 
dNf{b) 



PTdpTd(j) 



F,{k,b,(^)®{TS + SS) , 



(51) 



where the superscript sh denotes distributions involving 
shower partons. 



C. Momentum Shift 

Before proceeding to the calculation of {px , b) , let 
us first determine the normalization of the inclusive dis- 
tribution averaged over all 4>. It follows from Eqs. (j42p . 
and (HH) that 



diVf (b) 
PrdpT 



1 

2^ 



PTdpTd(j) 



—go{k,b)Mk)®iTS + SS) . 
Ztt 



(52) 



where the dependence in T is neglected, since this equa- 
tion will be used only for 0-10% centrality in this subsec- 
tion. For the TT component we take the average of B 
and B + R, and get 



dN^b) 
pxdp- 



cHb) 



n I 



}(53) 
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(GeV/c) 

FIG. 6: Fit of the tt'' distribution at one point, pr ~ 0.435 
GeV/c, by adjusting the energy loss coefficient eo- The data 
are for 0-10% centrality from Ref. [4^ . 

The sum of these two equations gives the total inclusive 
pion distribution dNT^{h)/pTdpT- It is to be compared 
with the data, labeled dN/2TrpTdpT, where N refers to 
the number of pions per event integrated over all (j). We 
have one parameter in Eq. (|52p . which is eg in Eq. (j37p . 
It determines the scale of momentum shift due to en- 
ergy loss of the hard parton in the medium, and enters 
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Eq. (|52|) through go{k,b). The larger eo is, the smaller 
is go{k, b) and the more suppressed is the shower parton 
contribution. We determine eq by fitting the normaliza- 
tion of the data for 0-10% centrality. Note that the pr 
dependence is not adjustable, but the overall degree of 
suppression of the shower component is adjustable, since 
eo is not determined within our formalism. Fitting the 
normalization at just one point (px = 4.35 GeV/c) we 
obtain 



0.55GeV^/2 . 



(54) 



The rest of the pr distribution is shown by the solid line 
in Fig. 6 in excellent agreement with the data [i^ for 
y/sNN = 200 GeV, the same set as used in . 
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FIG. 7: (Color online) Factional momentum shift for various 
normalized impact parameter b — b/2RA and (j). 

Since we have obtained fc'(fc, b, (p) in terms of t£ al- 
ready, we can determine the fractional shift Ak/k, shown 
in Fig. 7, for three values of normalized impact parame- 
ter, 6 = 0, 0.4, and 0.7, and three values of at 0, 7r/4, 
and 7r/2. They all decrease with k roughly as 1/Vk, as 
expected. There is, of course, no dependence on cj) when 
6 = 0. As 6 increases, the nuclear overlap region gets 
smaller, and Ak/k becomes smaller. At fixed 5 > 0, 
Ak/k is smaller at (/) = (in solid lines) than at </) = 7r/2 
(in dashed-dot lines) because the average path length is 
shorter. Since /3 in Eq. is around 8, G(fc, b, cj)) can 
be quite small, if Ak/k is not infinitesimal. We show 
in Fig. 8 the harmonic components gQ{k,b) and g2{k^b) 
defined in Eq. (gH) for 6 0, 0.4 and 0.7. The depen- 
dence of the hard parton is important because it gives the 
dominant contribution to the hadronic U2 . As 6 increases 
gQ{k, b) increases rapidly, indicating less suppression. For 
the lines showing g2{k,b), which are amplified by a fac- 
tor of 10 in Fig. 8, we see that their values are between 
0.03 and 0.04, relatively insensitive to fc, and, of course, 
(72(^,0) = 0. It is the ratio g2{k,b) /gQ{k,b) that sets the 
scale for V2 of hard partons. 




^ k (GeV/c) ^ 



FIG. 8: (Color online) The functions go{k,b) and g2{k,b) for 
some typical values of b. 



D. V2 for Pion 

We are now ready to calculate f J(pt, b). For the con- 
tribution from the shower partons at intermediate px we 
write out Eq. ([5T|) fully as 



dN': 



sh 



PTCtpTdcf) 



dkkF,{k,b, (f>) 



rPT 
XI dq 
/o 



PT -q 



(55) 



Putting this in Eq. jH) and using ([44]) and jHI) for 
the (p dependences, we perform the integrations over <^ 
first, and then over q and k. The result is fJ(pT,fe) for 
PT > 3 GeV/c. That is shown in Fig. 9 for the high 
Et portion in that figure. We note that the effect of (f> 
dependence of T'{b,(j)) on the magnitude of V2'^^{pT,b) 
in the high Et region is only about 10%. Thus the main 
contribution to V2'^^{pt, b) comes from the p dependence 



of Fi{k,b,(j)) for the hard parton in Eq. ([55]) . not from 
that of the thermal parton T in ([5T|) . If we neglected the 
4> dependence of T'{b, 4>) in ((55)) . we would have a simple 
formula in closed form 



",sh 



g2{k,b)h{k)®{TS + SS) 
go{k,b)Mk)^{TS + SS)' 



(56) 



which offers a succinct exhibit of all the important factors 
for the main contribution to V2 from the hard parton at 
high Et- 

For the low pT part the contribution from thermal par- 
tons is what we obtained in Sec. III. A, which we can write 
out explicitly from Eqs. Q and pi]) 



7r,th 



sin2$(5) 



tt/ [eET{PT)/n'-i^ + 2$ (6) 



(57) 



11 



0.3 



PHENIX 
40-60% 
5-1 0% 
0-5% 




2 3 
(GeV) 



FIG. 9: (Color online) for a wide range of Et- The small 
symbols are the same as those in Fi g. 1 ; the larger symbols 
(in green) are preliminary data from [47l |. 



The result has been shown in Fig. 1, and is shown again 
in Fig. 9 for Et < I GeV. 

The high and low Et regions should be connected 
by an interpolating function that spans the two regions. 
That function should depend on the relative weight of 
the pt distributions given in Eqs. ([52]) and ((53|) . For no- 
tational brevity let us use TT to denote dNf{h)/pTdpT, 
and TS + SS to denote dN^^{b)/pTdpT, which is the av- 
erage of Eq. ([55)1 over 2tt of (j). During the interpolation 
procedure that equation should not be used for pT too 
small. To ensure that, we cut it off by hand by inserting 
a cut-off function, 1 — exp(— 0.1p|n), which decreases from 
^ 1 to ^ as Pt is decreased from 3 to 1 GeV/c. This 
damping factor has been used in Fig. 6. We construct 
the weight function 



WipT,b) = 



TT 



TT + TS -t- SS 



(58) 



in terms of which we define the overall by 



{pT,b) = V2'^^{pT,b)W{pT,b) 



+vr\pT,b)[l-WipT,b)] . (59) 

Our calculated results for V2{pT,b) are shown in Fig. 9 
in terms of Et for < Et < 5 GeV at various values 
of b in fm, chosen to correspond to the centralities 0- 
5%, 5-10%, 10-20%, 20-30%, etc., as in Fig. 1, with c 
being approximately b^. It is evident that in all cases 
increases at low Et and decreases at high Et, reflecting 
the different mechanisms responsible for anisotropy. 

The data for Et < 1 GeV in small symbols are from 
Ref. already shown in pg. 1. For Et > I GeV we 
show the data from Ref. [47| for 3 centrality bins: 0- 
5%, 5-10%, 40-60%; they are indicated by large symbols. 
The first ones (0-5%) in large open circles correspond 
to b = 2.3 fm; the next (5-10%) in large filled circles 



correspond to & = 4.2 fm. The agreement with our results 
are good. There are no more centrality bins in Ref. [i^l 
that correspond to what we have calculated. We show 
only 40-60% in large triangles, which correspond to b — 
10.2-11.3 fm among our calculated curves. On the whole 
our result reproduces the data quite well. 



E. V2 for Proton 

For proton production we consider TTT, TTS and TSS 
recombination, leaving out SSS which is not important 
unless Pt > 8 GeV/c. The complication of 3-quark re- 
combination has been treated in Ref. [2^ already; it does 
not affect the calculation of V2, since hadronization oc- 
curs after jet quenching as we have seen in the pion case. 
Thus we proceed as before using Eq. dZl) for v;''^\pT,b) 
for Pt < 1 GeV/c with appropriate Bp{pT,b)/ Rp{pT,b) 
discussed in Sec. III.B. At intermediate pT we generalize 
Eq. (|51[) by including one more thermal parton, i.e.. 



dN^^ib) 
PTdpTdcj) 



F,{k,b, 4>)®{TTS + TSS) 



(60) 



When written out in full, it looks like Eq. ([55|) but with 
an extra factor, T(qi,6, 0), attached to each of the two 
terms in ([55]). For those thermal partons we use Eqs. 

and (gSl), except that T' in (gHl) has a b dependence, 
given m The calculation of wf' (pt, b) can then be 
carried out as before, but with an extra integration over 
qi . Although there is no analytical formula for the result 
of that calculation, a simplified form that neglects the 
dependence of T"(6, <j)) is 



{PT,b) 



g2{k,b)Mk)®iTTS + TSS) 
go{k,b)Mk)<E>{TTS + TSS) 



(61) 



which is the counter part of Eq. ([5S)) . The quantitative 
results shown below is, however, from the full calculation. 

The overall v^ijpT, b) for all pt region again involves a 
weighted average 

+vr''{pT,b)[l~W[pTM ■ (62) 
where the corresponding weight function is 

TTT 



W{pT,b) 



TTT + TTS + TSS 



(63) 



For TTT we have 



PTdpT nrriT 

while for TTS and TSS we have 



dK 



sh 



Prdpi 



sS V. 

[ dkkMk)go{k) 



X / dqidq2{TTS + TSS) 



(65) 
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where the average T over ( 

f'(6,0) = r'(&). 



is used, which means setting 



between our results and the i'2 data, on the whole we 
have been able to reproduce the data very well. 
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FIG. 10: (Color online) for a wide range of Et- The small 
symbols are the same as those in Fig. 4. The larger symbols 
(in green) are preliminary PHENIX data [47]] for centralities 
0-5% (open circle), 5-10% (fuU circle), and 40-60% (triangle). 
The STAR data are for 10-40% (blue square) at = 62.4 

GeV Si. 



The result on V2{ET,b) is shown in Fig. 10. The 
curves exhibit similar patterns as those for (-Et, b), but 
are generally higher, and persist to rise at high b when 
Et < 1.5 GeV. Although the (j) dependence of the hard 
parton is the same for tt or p, the thermal partons are 
not the same for the two cases. Since two thermal par- 
tons can participate in the formation of proton, the 
dependence of T'{b,(j)) plays a larger role in V2'^^{pT,b) 
than in v^'^^ipr, b). The data in the higher Eq^ region are 
shown by the same symbols as in Fig. 9. For 0-5 % and 
5-10% centralities the data are from Ref. 1471 for p + p. 
However, for 40-60% centrality the data in [43| for p -I- p 
differ significantly from the data in [4^ for p alone. It 
is the latter in large triangles that we show in Fig. 10. 
We include in that figure also the STAR data 0] at 62.4 
GeV for p + p a.t 10-40% centrality that agree well with 
our dotted line for 20-30%. For all the centralities shown 
the data are in general accord with our results. 

To summarize this long section we point out first that 
there is only one free parameter eg used to determine the 
medium suppression of hard parton by fitting the pion 
distribution at one point, pr = 4.35 GeV/c in Fig. 6. 
The parameter eo replaces the suppression factor ^ used 
earlier [2^ for central collision only without considering 
explicitly the issue of energy loss. The shape of the pr 
distribution in Fig. 6 is calculated. The path length de- 
pendence of jet quenching introduces the dependence 
at high pt that is the main source of elliptic flow. The 
characteristics of U2(£'t,&) for h — n and p in Figs. 9 
and 10 are determined without any more adjustable pa- 
rameter. Although there are some minor discrepancies 



V. BREAKING OF QUARK NUMBER 
SCALING 

In the naive application of the recombination model 
there is quark number scaling (QNS) of V2, which may 
be stated as the universality of V2{pt /n'q)/ng where n„ is 
the number of constituent quarks in the hadron h [20|.|45|. 
The simple argument used is based on the factorizabil- 
ity of the distribution of the quarks that recombine, and 
on the simplification of the recombination function to the 
form that contains d{qj —pr/nq) for each of those quarks. 
There are also other considerations for the origin of QNS 
(49| . Experimental verification of QNS has evolved to 
the replacement of pT by Er with impressive confirma- 
tion of the scaling behavior [1]- [III : at least at low Et / Uq . 
Since it is known that fragmentation is more important 
than recombination at very high pT (or in very periph- 
eral collisions), QNS should break down at some point 
[ll,[l^. The question is at what point. We show here 
that it occurs rather early, even when TS recombination 
is still dominant. In fact, at even lower pT where TT 
and TTT recombination are more important, QNS is not 
valid in general for specific centralities, as discussed in 
Sec. III.B. However, averaging over all centralities leads 
to approximate QNS, in agreement with minimum bias 
data [a, J3.,.ll, 12). 

The breaking of QNS is due mainly to the breaking of 
factorization of joint parton distribution, which, if true, 
would have the form 



En, • • • ='[[F^{q^,<j>i) 

1=1 

rig 

= Y[F,iq,){l + 2vliq,)cos2^,) . 
i=l 

From this follows 



V2{pt) - nqvl{pT/nq) 



(66) 



(67) 



for qi — pT/nq. For pion at low pt the factorization 
of qq distribution is reasonable, and that is how Eqs. 
([TT|) and (|12p are obtained at all centralities. For pro- 
ton at low Pt the factorization of uud distribution be- 
comes questionable for non-central and peripheral colli- 
sions, the consequence of which is that the inverse slope 
for proton differs from that of the single inclusive u (or d) 
quark. Instead of investigating the non-factoriable form 
of Fuudiqi, q2, qa), we have adopted in Sec. III.B phe- 
nomenological form for the average of Bp{pT) + Rp{pT , 4') 
over all </> with given in Eq. ([^5]) and shown in Fig. 3. 
That results in the centrality dependence of given in 
Eq. (PT|) . Since T^' is constant for all c at the value speci- 
fied in Eq. V2{pt-, b), as given by Eq. ([5]), has initial 
Et dependences that are different for pion and proton. 
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as expressed explicitly in Eqs. ((22)) and (|32p . except at 
one point of cross-over that occurs at a centrality roughly 
equivalent to minimum bias. 

At intermediate px where the contribution from shower 
partons is important, the (jj dependence comes mainly 
from Fi{k,b,(f>) for the hard parton in Eq. ([5T|l for 
pion and Eq. (|60p for proton. Thus the corresponding 
V2{pT,b) for h — TT and p are given primarily by Eqs. 
(|56)) and (|6T|) that have the same structure, the com- 
mon factors of importance being 172(^1 b) in the numera- 
tors and go{k, b) in the denominators. The significance of 
g2{k,b)/ go{k,b) being the major factor for both V2 and 
V2 immediately implies that 



^vUET/2,b)>^v^2(ET/3,b) 



(68) 



at intermediate Ex for all b. However, when the (f> depen- 
dence in T in those two equations is taken into account, 
the results on v!^'^^ are enhanced more for proton than 



for pion. That makes the inequality in Eq. (j68p less un- 
balanced. Nevertheless, QNS is broken because not all 
partons contribute equally to the azimuthal anisotropy, 
even though both hadrons are formed by recombination. 

Since Figs. 9 and 10 exhibit all properties of t;J(-ET, b) 
and V2{ET,b) that we have obtained, we can calculate 
v^{ET/2,b)/2 and v^{ET/3,b)/3 and show their differ- 
ences for four centralities in Fig. 11. Although the pion 
(dashed) and proton (dashed-dotted) lines start out with 
similar slopes at low Et, they separate at higher Et, the 
former being consistently higher than the latter. The 
data that support QNS are for minimum bias and at 
low Et where the theoretical lines are in agreement with 
the data at all four centralities. The breaking of QNS 
becomes visibly clear at Er/nq > 0.5 GeV, and is dis- 
cernible even in the minimum bias data of Ref. Q . 

The reasons for plotting V2{Et /nq)/nq are not only 
to check QNS, but also to give a hint of the elliptic flow 
of the recombining quarks, if QNS were valid. We can 
calculate u^C-E't, b) for quarks theoretically before recom- 
bination. At low Et Eq. ([7]) is valid for quarks also, with 
Bq{q)/Rq{q) being given by Eq. (|20)). For Et we use the 
constituent quark mass, which we take to be 0.3 GeV. 
The resulting V2{ET,b) may be regarded as reliable for 
Et < 0.7 GeV. For higher Et we recall Eqs. ([51]) and 



(|56p for pions. The structure for quarks is similar, ex- 
cept that there is no recombination. Thus for a quark of 
type j in the shower we have 



dNf{b) 
qdqdcj) 



j dkkF,{kA <t>)Sl{q/k) 



^,,sh. /rffcfc.g2(fc,b)E./.(fc)^f(g/fc) 

S dkkg^{k,b)Y.J,{k)Sl{q/k) 



(69) 



(70) 



These equations are for shower partons with momentum 
q not too low. The lower bound for their validity may 
be set at q ^ 1.5-2.0 GeV/c. For illustrative purpose 
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FIG. 11: Scaled 112 for pion, proton and u quark. For the u 
quark the thermal «2 is plotted at low Et, changing to shower 
V2 at high Et- The data points are scaled from those in Fig. 
1 (open circles) and Fig. 4 (filled circles) from Q. At higher 
Et /uq only minimum bias data are available, not suitable for 
display here. 



we consider u quark for definiteness. For each centrality 
in Fig. 11 we show W2 ((?,&) for the two regions: thermal 
quark for Et < 0.7 and shower quark for Et > 1.5 GeV; 
in between we simply connect the two by a smooth inter- 
polating curve without mathematical significance. The 
physical significance of V2{ET,b) shown is that it is aris- 
ing function of Et at low Et and a gently falling function 
at higher Et, characterizing the thermal and shower par- 
tons, respectively. 

It is now evident from Fig. 11 that the three curves 
for u, TT and p for each centrality are distinctly different. 
From central to mid-central collisions up to 40% central- 
ity the difference between tt and p curves may not be 
large enough for the present data on V2 to discriminate. 
There is some evidence for QNS breaking at Et/uq > 0.5 
GeV in the minimum bias data for tt and p [1J|. They 
are for ^/snn = 62.4 GeV. Roughly, the pion data are 
above the proton data, as in Fig. 11. Current experimen- 
tal efforts have largely been to make corrections due to 
fluctuations in eccentricity in order to achieve QNS. Our 
suggestion is to focus on individual centrality bins and 
quantify the breaking of QNS. 

The results shown in Fig. 11 are from extensive and 
detailed calculation. It would be helpful if the quantita- 
tive difference between the two behaviors for pion and 
proton can be explained in simple terms, albeit inexact 
and schematic. To that end let us assume dominance 
of TS and TTS recombination for ^'r/n, > 1 GeV and 
write 



(cos 20) 



TS ' 



(cos 2c, 



'TTS 



(71) 



If the contribution to V2 from thermal and shower partons 
are expressed as vj and ff, respectively, then we have 



approximately 

v^~vj + vl , vl ~ 2vJ + vl , (72) 

where the moraenta of the quarks are not shown exphc- 
itly because they are comphcated to express. Unhke the 
naive description given in Eq. (67), the thermal parton 
momentum is smaller compared to the shower parton mo- 
mentum because the distribution of the former is damped 
exponentially while that of the latter is power suppressed. 
Furthermore, even if the quark momentum is q — px/nq, 
we have vj (q) < i;| {q) because of the difference in the (j) 
dependence in T and S. That inequality becomes more 
unequal when smaller q is used in vj{q) and larger q in 
V2 (q) ■ Thus roughly we have at some fixed Et > 2 GeV 

vUEt) ^ 2 + J ^ 2 

where 

6 = vliq+)/vJ{q-) - 1 > 0, q±> Et/2 . (74) 

It is now clear that at the very basic level the breaking of 
QNS is due to the non-equivelance of the (p dependences 
of thermal and shower partons. Experimental verification 
of the QNS breaking can therefore render an indirect sup- 
port for the recombination mechanism involving thermal 
and shower partons, as we have described. 

VI. CONCLUSION 

At y/SNN = 200 GeV the density of partons with mo- 
mentum fraction x ^ 0.03 is high, and their scattering 
into fcT 3 GeV/c can readily occur in a heavy-ion colli- 
sion. Thus the formation of ridges due to weak jets is an 
aspect of the event structure that is pervasive, and should 
be taken into account in the study of distribution. Hy- 
drodynamical expansion of the dense medium undoubt- 
edly takes place, but there is no requirement from any 
fundamental principle that the thermodynamical descrip- 
tion must be valid within 1 fm/c after collision. With 
semi-hard scattering driving the azimuthal anisotropy, 
fast equilibration need not take place, but hydrodynam- 
ics can still have its realistic application at later time. 
Short time-scale physics implies hard or semi-hard pro- 
cesses by uncertainty principle, and it is not the proper 
domain of equilibrium physics. We have demonstrated in 
this study that the observed features of elliptic flow are 
on the whole reproduced by ridge consideration at low pT 
and thermal-shower recombination at intermediate pr- 
The validity of hydrodynamics at late time is implicitly 
incorporated in our approach when the bulk parton dis- 
tribution is taken to be thermal, i.e., exponential before 
hadronization. How to combine semi-hard scattering and 
hydrodynamical expansion is a time-evolution problem 
worthy of careful investigation. 

Although we have used recombination in all px regions, 
we have shown that quark number scaling is not generally 
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valid. At low Et it is approximately valid because the 
partons that recombine are mostly multiplicative, each 
having roughly the same V2- At intermediate Et the 
thermal and shower partons contribute to V2 differently, 
so the mechanism that gives rise to QNS is lost. 

Since rapid thermalization is not required in this study, 
the partons in the dense medium need not interact more 
strongly than in the usual formulation of QCD. Viscosity 
need not be low, since hydrodynamical calculation should 
be redone with different initial conditions. No where in 
this study suggests that the medium created in heavy-ion 
collisions behaves as perfect fluid. 
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APPENDIX A: ANGULAR DISTRIBUTION OF 
RIDGE PARTICLES 

We describe in this Appendix the subprocesses in- 
volved that lead to the distribution of particles in the 
ridges. The subject of correlation that depends on the 
trigger direction is outside the scope of the present prob- 
lem, and will be described separately in [35|. Our in- 
terest here on the single-particle inclusive distribution is 
less complicated; however, it does require an integration 
over all angles of the semi-hard partons that generate the 
ridges. Thus some aspect of the correlation problem will 
be adopted here from (35| . Since the problem of angular 
correlation does not depend critically on the magnitude 
of pt, we suppress the pt variable in this Appendix until 
the very end where contact with the main text is to be 
made. 

It is important to recognize two time scales in this 
problem: one is the semi-hard scattering at early time 
that is sensitive to the initial configuration of the system, 
and the other is the hadronization process at later time 
for which the elliptic geometry is more relevant. For the 
former the almond-shaped region is bounded by two cir- 
cular arcs of radius Ra with centers at a; = ±6/2, y — 0, 
where the (x, y) coordinates are centered at the origin of 
the almond with the x-axis on the short side and in the 
reaction plane. We refer to it as the A geometry. In this 
Appendix we normalize all lengths by RA^ so the radii 
are 1. For the latter the ellipse is defined by 




where w = 1-5/2, h = [l-(6/2)2]i/2^ and u = 1 at initial 
time. We refer to this as the E geometry. We assume 
that, as the dense system expands, we need only let u 
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increase in using (|Aip to describe the boundary of the 
system. The vector normal to the surface at any point 
{x,y) is the gradient of u, whose azimuthal angle is 



'ip{x,y) ~ tan 



\h) X 



(A2) 



It specifies the direction of local flow. 

Semi-hard scattering can occur at any point inside the 
almond. Let it be labeled by P whose coordinates in the 
A geometry are {x,y). Limiting ourselves to only the 
2D transverse plane on the basis that all particles de- 
tected in the final state are in a narrow rapidity bin at 

= 0, we consider one of the created semi-hard partons 
moving in a direction at azimuthal angle (ps toward the 
boundary, away from the interior. Let t be the distance 
between P and the boundary measured along the trajec- 
tory of the parton. Because of energy loss to the medium, 
t cannot be too large if the parton is to emerge and give 
rise to a ridge in addition to a trigger particle. In |35j] 
a distribution in t is considered with parameters deter- 
mined by phenomenology. For our purpose here that is 
more illustrative than data fitting, we adopt the discrete 
average value of t = 0.1. That is, we assume that all 
semi-hard partons are created along the inner circles of 
radius = 0.9 from either center in the A geometry. Let 
us refer to the locus of those points as the rim. Thus P 
is a point on the rim, and can be specified by the angle 
9 measured from the pertinent center. For definiteness, 
we consider the right half of the almond rim, so 9 varies 
from — <I>o to -t-'&o, where $0 = cos~^(6/2ro). Recoil par- 
tons that move toward the interior, as well as those that 
are created there, are totally absorbed and contribute to 
the bulk. High momentum jets can get out in any direc- 
tion, but they are rare and do not contribute to the ridge 
particles that can influence the azimuthal anisotropy. 

The process of energy loss to the medium by the semi- 
hard parton cannot be calculated reliably. In 35] suc- 
cessive soft emission with probability proportional to the 
local density is modeled with correlation to the jet trajec- 
tory. Conversion of the lost energy to the enhancement of 
thermal partons, which subsequently hadronize to form 
ridge particles, is also hard to treat rigorously. It is, 
however, reasonable to assume that for every semi-hard 
parton at (ps there is a corresponding cluster of ridge 
particles in </> that has a Gaussian distribution around 
(ps- This is essentially a description of what has been 
observed in ridge phenomenology (23 |. where the trigger 
direction takes the place of 4>s here. More specifically, we 
consider the correlation function 



where 



C{x, y, (j), cf)s) = D{x, y)G{(j), (f>s 



G{cb,(bs)=exp [-(</) -0,)V2A] 



Dix.y) = TaIs) 



1 -e 



^crTB{\s-b\) 



(A3) 



(A4) 



(A5) 



+ TBi\s-b\) [l-e-"^^(^) 



The last quantity D{x,y) is the local nuclear density at 
(x, y), related in the Glauber model to the thickness func- 
tion Ta{s) and pp inelastic cross section a in the conven- 
tional way that will not be detailed here. We are not con- 
cerned with the overall normalization because the aim of 
this Appendix is to derive the (p dependence of the col- 
lected ridge particles generated by semi-hard partons at 
all points along the rim, and the normalizations of the 
ridge versus bulk are considered along with the pr de- 
pendence in Sec. III. The parameter A in (|A4p is found 
in [3^ to be 0.11, which we use in the following. It cor- 
responds to a half- width of about 20°. 

The semi-hard parton is scattered to an angle cps ran- 
domly, so the inclusive distribution of the ridge particles 
without trigger must average over all cps . This should be 
done for each point P on the rim. A mathematically sim- 
ple and physically reasonable approximation of the result 
of that averaging is that the average direction of the ridge 
particle, 4>, is normal to the surface, since that is the only 
direction in the problem. That is the approximation used 
in [23| . where only the A geometry is considered. Here 
we treat the problem more quantitatively by averaging 
over <ps at each (x, y) point and then integrating over all 
points along the rim. Furthermore, we improve the cal- 
culation by using E geometry, since hadrons are formed 
later when the system develops elliptical shape and they 
follow the direction of the local flow. Although G(0, (ps) 
in (jA4p makes no reference to the system shape, the av- 
eraging process is sensitive to it. That is shown explicitly 
as follows 

dN 1 

-n:i^'y)^o d(l)sC{x,y,<p,(ps) , (A6) 

where the range of integration is ±1 around the normal to 
the elliptic surface ip{x,y), given in (|A2p . which specifies 
the direction of flow at {x, y) near that surface. We have 
considered wider range, but the difference from the above 
is negligible. The averaging in (jA6| accounts for the semi- 
hard partons that move toward the surface and emerge 
from it. The difference between the surfaces in the A 
and E geometries becomes significant, when P is near 
the tips of the almond region. However, the density is 
lower there, so the effect of the difference on the yield is 
not prominent. 

To see the results at various points on the rim, we show 
in Fig. 12(a) the azimuthal distribution of the ridge parti- 
cles for b = 0.6 created by semi-hard partons, originated 
at three illustrative points, whose azimuthal angles, mea- 
sured from the center of the right half of the rim, are 
denoted by 9, where 9 = tan"i[y/(a; -I- 6/2)] . The thin 
solid line exhibits the shape of G{(f>, (ps) for (ps = 0. The 
thicker solid line shows dN{x, y) / d(p for = 0. Evidently, 
the averaging process widens the <p distribution, since jets 
at any angle in the range Q < \(ps\ < ^ can contribute. 
The other two curves in that figure show the contribution 
from points at 6* = 25° (dashed) and 50° (dash-dotted). 
Their peaks are lower because of the D(x,y) factor in 
[3I) and are centered around <p = 9. The experimental 
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FIG. 12: (Color online) (a) Azimuthal distributions of ridge 
particles due to semi-hard partons created at three points on 
the rim mt 9 ^ 0° (solid), 9 = 25° (dashed), and 9 = 50° 
(dash-dotted) for b = 0.67?a- Thin solid line is the Gaussian 
correlation function G{<j),4's) for <l)s =0°. (b) Azimuthal dis- 
tributions of ridge particles after integration over all points 
on the rim. b is labeled in unit of Ra- 



detector cannot distinguish the different points of semi- 
hard scattering, so we must integrate over all 6 



(A7) 



The result for 6 = 0.6 is shown in Fig. 12(b). We note 
that the long tail on either side extends well beyond cj) = 
±90° . It is a consequence of the E geometry that we have 
used, since the normal to the ellipse near the top can be 
^p « 7r/2. Thus at the tip of the rim hadrons formed 
from semi-hard partons on the two sides of the rim can 
overlap, resulting in doubling the value of dN/d(l)\^=Tr/2 
that is calculated from (jA7[) . which integrates only the 
contribution from the right side of the rim. 

In Fig. 13 we show the 4> distributions for four values of 
b in the range — 7r/2 < (j) < n 12. In each panel the main 
contribution (in dashed line) comes from the right side of 
the rim where B is between — $o(^) and -t-$o(^) as seen 
in (jA7p . The dashed-dotted lines show the contributions 
from the left side of the rim — from Q > vr — $o(&) as 
well from 6 < tt + ^o{b) measured from the center of the 
left half of the rim. The sum is shown in solid (black) 
line. We have checked that at & = the distribution is 
absolutely flat, as it should be. The normalization will 
be discussed presently. 

For the purpose of rendering simple mathematics and 
transparent physics for the determination of V2{pt) in 
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FIG. 13: (Color online) Azimuthal dependences of ridge par- 
ticles for b = 0.3, 0.6, 0.9 andl.2 in units of Ra- Solid (black) 
line is the sum of the contributions from the right-side rim 
(dashed line) and the left-side rim (dash-dotted line). Straight 
(red) line is the box approximation with width at \<j)\ = ^{b). 
Vertical scale is adjusted such that the box height for b = 0.6 
is 1; scales for all other panels are not readjusted. 



Sees. II and III, we now approximate ^{b) shown in 
Fig. 13 by step function shown in solid straight (red) line, 
whose width is given by Eqs. (1) and (2). The values of 
$(6) for b = 0.3, 0.6, 0.9, and 1.2 are 81.4°, 72.5°, 63.3° 
and 53.1°, respectively. The height of each inverted box 
is determined by matching the area under the box with 
the area under the calculated ridge curve, which follows 
from (A3) through (A7). We have given those equations 
arbitrary overall normalization because the true normal- 
ization of the ridge distribution, R(pT,4'), involves the 
consideration of pT dependence, which is discussed in 
Sec. III. In Fig. 13 we have adjusted the vertical scale so 
that the height of the box approximation for b = 0.6 is 1. 
The relative heights among the different panels of that 
figure are not adjustable. Thus Fig. 13 exhibits clearly 
the comparison between the results of the calculation of 
the (j) dependence and the box approximations for various 
values of 6, showing that the heights of the boxes remain 
essentially the same and that the widths specified by <&(&) 
summarize the effective dependence on centrality. In the 
lower-right panel of Fig. 13 for 6 = 1.2 the two wings of 
the yield curve extend considerably beyond the box be- 
cause of the difference between the A and E geometries 
when (j) is near 7r/2. The rapid descends of the wings are 
acceptably represented by the narrowing of the width of 
the box. 

Although details in each subprocess considered in this 
Appendix can be made more elaborate with more pa- 
rameters, the general property of the outcome is clear 
and cannot deviate too much from what is shown in Fig. 
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13. The box approximation given in Eq. (4) captures the and faciUtates the derivation of the very simple, analytic 
essence of the (j) dependence of the ridge yield, R{pt, <I)), formula, Eq. (8), for V2{pt, 
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